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Abstract 



We introduce a model of a two-dimensional (2D) optical waveguide with 
Kerr nonlinearity, linear and quintic losses, cubic gain, and temporal-domain 
filtering. In the general case, temporal dispersion is also included, although 
it is not necessary. The model provides for description of a nonlinear planar 
waveguide incorporated into a closed optical cavity. It takes the form of a 2D 
cubic-quintic Ginzburg-Landau equation with an anisotropy of a novel type: the 
equation is diffractivc in one direction, and diffusive in the other. By means of 
systematic simulations, wc demonstrate that the model gives rise to stable fully 
localized 2D pulses, which arc spatiotemporal "light bullets" , existing due to 
the simultaneous balances between diffraction, dispersion, and Kerr nonlinear- 
ity, and between linear and quintic losses and cubic gain. A stability region of 
the 2D pulses is identified in the system's parameter space. Besides that, we 
also find that the model generates ID patterns in the form of simple localized 
stripes, which may be stable, or may exhibit an instability transforming them 
into oblique stripes with zigzags. The straight and oblique stripes may stably 
coexist with the 2D pulse, but not with each other. 

Keywords: light bullet, zigzag pattern, diffraction, dispersion, filtering, non- 
linear amplification 
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1 Introduction and formulation of the model 



Equations of the Ginzburg-Landau (GL) type constitute a class of universal 
models to describe pattern formation and spatiotemporal chaos in nonlinear 
media combining dissipative and dispersive properties 0. The GL model of 
the cubic-quintic (CQ) type was first introduced by Petviashvili and Sergeev 
[H in the two-dimensional (2D) case, with an objective to construct stable fully 
localized 2D objects. In several works, stable localized pulses have been obtained 
in the cubic-quintic GL equation in one dimension (ID) [||. 2D pulses in models 
of the GL type were considered by Thual and Fauve and by Deissler and Brand 
[Q. Further, Firth and Scroggie predicted localized pulses in a two-dimensional 
GL equation with the nonlinearity corresponding to a saturable absorber in an 
optical cavity ||^. Localized 2D pulse was also found in the CQ equations of 
the Swift-Hohenberg type . Recently, a stable localized spiral pulse has been 
found in the 2D GL equation with the CQ nonlinearity 0. 

In previous works (^-0, localized pulses were sought for in isotropic 2D 
equations of the GL type. In this work, we aim to introduce a specific 2D 
GL equation which governs spatiotemporal evolution of light in a planar active 
(i.e., equipped with an intrinsic gain) nonlinear optical waveguide, and find sta- 
ble fully localized pulse solutions to it. The spatiotemporal 2D GL equation for 
the optical system is naturally anisotropic, as the temporal variable plays the 
role of one of the coordinates, see below. The existence of stable fully localized 
pulse solutions to this equation suggest a new experimental realization of "light 
bullets" , i.e., spatiotemporal optical pulses. Thus far, "bullets" were experimen- 
tally observed only in conservative second-harmonic-generating optical media in 
an (effectively) 2D geometry We aim to propose a medium in which truly 
robust anisotropic 2D "bullets" exist as a result of stable balance not only be- 
tween the Kerr nonlinearity, spatial diffraction, and temporal dispersion, but 
also between various loss and gain mechanisms. 

A prototype model equation, which governs the evolution of a local ampli- 
tude u(z,x,t) of the electromagnetic field in a 2D waveguide, is 

iuz + ^u^x + ^ (/3 - iF) Utt + lupM = ijou. (1) 

Here, z and x are the propagation and transverse coordinates, 

T = t- z/Vo (2) 

is the so-called reduced time, where t is the physical time, and Vq is the group 
velocity of the carrier wave. The cubic term and the one ^ u^x m Eq. (^ 
account for the self-focusing Kerr nonlinearity and transverse diffraction of the 
light pulse, respectively. The linear terms with the coefficients /3 and F take into 
regard chromatic dispersion and optical filtering, which is formally tantamount 
to diffusion in the r-direction. The term on r.h.s. of Eq. (|l|) with 70 > 
is the optical gain (amplification). The gain is necessary to compensate losses 
introduced by the filtering. Note that filtering may be induced either by optical 
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filters directly inserted into the waveguide, or by a choice of an amplification 
mechanism which provides for a bandwidth-limited gain. 

Due to its meaning, the filtering coefficient F in Eq. is always positive, 
while the dispersion coefficient (3 may be negative (corresponding to normal 
dispersion), positive (anomalous dispersion), or (if the carrier wave is launched 
at the zero-dispersion point |^). In fact, just the latter case, /3 = 0, is the most 
fundamental one, as it corresponds to the simplest model of the present type 
that can give rise to 2D pulses, see below. 

The diffraction coefficient (the one in front of the term Uxx in Eq. (0)) in 
optical media has no imaginary part, hence Eq. (^ is intrinsically anisotropic: it 
features only diffraction in the x-direction, and, in the above-mentioned simplest 
case (3 — Q, only an effective diffusion (alias filtering) in the T-direction, while 
z plays the role of the evolutional variable. Other 2D anisotropic models of 
the GL type are known in superconductivity (see, e.g., Refs. 0), as well 
as in the general pattern-formation theory ||T^ - However, in models of 

superconducting layers anisotropy is very different from that in the present 
model; it is usually represented by an asymmetric diffusion operator in the GL 
equation proper, which is coupled to an equation for the magnetic field. In the 
above-mentioned pattern-forming models, anisotropy is also drastically different 
from what we have in Eq. (|l|); in those models, the anisotropy is induced through 
a relative asymmetry of the diffusion operators in two coupled GL equations 
[ p!2| , or by addition of extra terms which directly break the axial symmetry 
of the equation iQ, or if a higher-order equation of the Segel - Whitehead - 
Newell type is considered To the best of our knowledge, no 2D model in 
which the GL equation is diffractive in one direction and diffusive (or mixed 
diffusive/diffractive) in the other has been introduced before. 

Stationary solutions to Eq. (^ are sought for as 

u{z,x,T) = e'^''U{x,T), (3) 
where the complex function U{x,t) obeys an equation 

-kU+ hjxx + ^F) Urr + \U\^U = i-foU. (4) 

In particular, ID pulse solutions of Eq. (|l|) have the form U{x, t) — U{9), where 

6 = T ~ px, (5) 
with an arbitrary real constant p. In this case, the function U obeys an equation 

- fcf/ + i [/ + (/3 - tF)] ^ + \U\^U = ijoU, (6) 

and is assumed to vanish at — > ±oo. 

Equation (^ coincides with the stationary version of the usual ID complex 
cubic GL equation, and at all values of p it has a single solitary-pulse solution 
of the form 

C/sp(r)-^ [sech(r;0)]l+^^ (7) 
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where 



if = 27o[F-2(/3 + p2)^; 



I -1 



= i{ri/2Ff[{p + p')+F^] v/9C^ + p2)T4F2 - 3 (/3 + 

Thus, Eq. has the solution in the form of a ID stripe oriented in any direction 
in the (x, t) plane, except for that corresponding to p = 0, i.e., parallel to the r 
axis, see Eq. (||). This suggests to look for a solution in the form of a 2D pulse 
localized both in x and in r, i.e., satisfying the boundary conditions 

= oo) = [/(|t| = oo) = 0. (8) 

In the context of the 2D nonlinear optical waveguide, the 2D pulse ("light 
bullet") is traveling at the carrier- wave's group velocity Vq in the laboratory 
reference frame, see Eq. (||), therefore the most physically relevant case, in 
which a stable "bullet" may be observable, is that when it is traveling in a closed 
ring. The latter situation may be realized if the nonlinear planar waveguide is 
a part of a closed optical cavity (resonator) equipped with intrinsic gain (see, 
e.g., Refs. Il^, |l^), so that Eq. (|^) is an average equation for the cavity. A 
more exotic, but also more straightforward, possibility to realize the closed ring 
is to roll the planar waveguide described by Eq. (|l|) into a cylindrical surface, 
so that the propagation coordinate z will become a cyclic one. 

However, both ID and 2D pulses, that may be found as solutions of Eq. (^, 
are unstable as a solution to Eq. (|l|) , because the presence of the linear gain in 
this equation makes the zero solution, which is the background of the localized 
pulse, unstable. Thus, to obtain physically meaningful pulses, it is necessary 
to modify the model so that to make its zero solution stable. One possibil- 
ity is to linearly couple the cubic GL equation to an extra linear dissipative 
equation, following the pattern of similar ID models Physically, the 

accordingly modified model corresponds to a dual-core waveguide, in which one 
(nonlinear) core is active (has the gain), while the other one features only losses, 
the filtering and nonlinearity in it being negligible. Within the framework of 
the ID geometry, it has been demonstrated that, in a vast parameter region of 
the dual-core system, the zero solution is stable and, simultaneously, there are 
two stationary solitary pulses (they were found in an exact analytical form ||l9| , 
following the pattern of the exact solution (^). As a result, the pulse with the 
larger amplitude is stable |l^ . It is quite feasible that stable 2D pulses may 
also exist in a 2D version of the dual-core model. 

Another possibility to produce stable pulses is to insert directly into the 
waveguide, or into the above-mentioned closed cavity incorporating the waveg- 
uide, a nonlinear saturable absorber (NSA); note that optical resonators includ- 
ing gain (usually provided by the Erbium dopant), nonlinear waveguide, and 
NSA are well known (in the ID case) in the context of soliton-generating lasers, 
see Refs. |l^] and references therein. In the presence of NSA, it is easy to 
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select the system's parameters so that, for a weak light signal, the absorption 
is stronger than the linear gain, hence the zero solution is stable. However, as 
the absorption saturates with the increase of the signal's intensity, the system 
may feature nonlinear gain, opening way to the coexistence of the stable zero 
solution with stable pulses. 

As it is well known, the simplest model describing a combination of the linear 
gain and NSA is based on the above-mentioned GL equation of the cubic-quintic 
(CQ) type. In particular, the accordingly modified 2D equation (|l]) takes the 
form 

iUz + ^Uxx + ^ (/3 - iF) Urr + = -«7ou + «7i|MpM - i^2\u\'^u, (9) 

where the linear and quintic dissipation coefficients 70 and 72, as well as the 
cubic-gain coefficient 71, are positive. One may fix 70 = 1 and = 1 by means 
of obvious rescalings, then the eventual form of Eq. (H) is 

iUz + -Uxx + 2 (/5 ~ = - [lu + [l - 171) \u\^u + i'y2\u\'^u\ . (10) 

Strictly speaking, in the case P — (at the zero-dispersion point), which 
will be considered, among other cases, below, one should add the third-order 
dispersion (TOD) to Eq. (^, so that it will take the form §] 

iUz + -Uxx + -i) Utt ~ jiPuTTT ^ -[iu+ {l~ i^ii) \u\^u + i72|u|^u] , 
2 2 

where /3 is a real TOD coefficient. However, TOD is important only for pulses 
which are extremely narrow in the temporal direction, therefore this extra term 
will not be taken into regard in most cases in this work (but, nevertheless, the 
effect of this term on the 2D pulses will be considered below). The model may 
also include a quintic correction to the Kerr nonlinearity, but such a term it is 
not considered in this work. 

In the rest of the paper, an objective is to find, by means of direct simulations 
of Eq. (|l^), various stable stationary states in the present model, both 2D 
solitary pulses (corresponding to the spatiotemporaindocument "light bullets" 
in the waveguide model formulated above) and (quasi) ID states. The latter 
ones will be represented by stripe patterns which are localized in one direction, 
and by more sophisticated patterns that include steps (they are similar to zigzag 
patterns known in some other models I^^sults for the 2D pulses are 

presented in section 2, and ID patterns are described in a brief form in section 
3. Section 4 concludes the paper. 

2 Stable two-dimensional localized pulses 

Equation ( [lo|) was solved numerically by means of a pseudospectral code with 
256 X 256 modes and periodic boundary conditions in x and t, keeping fixed 
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sizes of the system in both directions, — Lr = L ~ 20. To generate the first 
stable 2D pulse, an initial configuration was taken in the form of a Gaussian, u 



2 exp 



— (^{x — + {t ^ /4 (i.e., the center of the pulse is placed 

at a; = r = L/2). In order to generate stability diagrams (see below), we 
then varied the parameters by small steps, the initial configuration for each 
simulation being the stable pulse produced by the simulation at the previous 
step. 



Figure 1: A typical example of a stable stationary two-dimensional pulse found 
in Eq. ^) at /? = 0, 71 = 2.5, 72 = 1/2. The absolute value of the field is 
shown as a function of x and r. 



A typical example of a stable fully localized pulse produced by the simula- 
tions is displayed in Fig. 1, which pertains to the simplest and most fundamental 
case, /3 = 0. To further illustrate the structure of the pulse, in Fig. 2 we display 
its field, represented in the form u{z,x,t) = a{x,T) exp[ikz + i4>{x,T)] with 
real amplitude a and phase cf), along two principal cross sections, x = L/2 and 
T = L/2 . It is noteworthy that the pulse is much narrower in the x (spatial) 
direction than in the temporal r-direction, which may be explained as follows: 
the diffusion (optical filtering) in the r-direction, which is a dissipative feature, 
produces an essentially stronger effect than diffraction (a conservative feature) 
along the x-direction; hence, to compensate the nonlinear terms, a much larger 
curvature of the soliton's profile at its center is necessary in the x-direction than 
in the direction of r. It is noteworthy that the so-called chirp, i.e., curvature 
of the phase profiles, which plays an important role in dynamics of solitons in 
nonlinear optics is strongly concentrated at the center of the pulse; note also 
a peculiar shape of the phase field in the x-direction. 

Results of many runs of the numerical simulations at different values of the 
control parameters were collected in the form of stability diagrams displayed in 
Fig. 3. They show regions in the parameter plane (71, 72), at fixed P — 0, and in 
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X. r X, r 

Figure 2: The amplitude (a) and phase (b) profiles of the same stable soliton 
which was shown in Fig. 1, taken along the cross sections t = L/2 (continuous 
lines) and x = L/2 (dashed lines). 



the plane (/3,7i), at fixed 72 = 1/2, in which stable 2D pulses have been found. 
The fact that the stability region expands while the dispersion coefficient 13 is 
changing from negative to positive, see Fig. 3(b), can be easily understood: the 
existence of a solitary pulse as a result of the balance between the self-focusing 
Kerr nonlinearity and anomalous dispersion (corresponding to /3 > 0) is much 
more natural than in the case when the self-focusing nonlinearity competes with 
normal dispersion (which corresponds to /3 < 0) [||. Nevertheless, stable pulses 
are found in the normal-dispersion area too. The pulse can exist in this case 
because the effect of the diffusion (filtering) in the r-direction beats the anti- 
localization effect exerted by the normal dispersion. 

To illustrate the change of the pulse's shape while varying parameters inside 
the stability region, in Fig. 4 we display its widths Wt and Wx of the pulse 
in the r- and a;-directions as functions of 71. In fact, the results presented 
in Fig. 4 pertain to a cut across the stability area delineated by Figs. 3(a) 
and 3(b). The cut is taken in the direction parallel to the 71 axis, at fixed 
values 72 = 1/2 and /3 = 0. The widths Wt and Wx were measured along 
the cross sections x = L/2 and r = L/2 of the pulse, as separation between 
points where, respectively, \u{x = L/2,t = L/2 ± Wr/2)\ — (l/2)|u|max, and 
\u{x = L/2 ± Wx/2,T = L/2)\ — (l/2)|u|i„ax, with |M|max taken at the center 
of the pulse, x = t = L/2. As well as in the case shown in Figs. 1 and 2, 
one concludes, looking at Fig. 4, that the pulse is essentially narrower in the 
x-direction. 

The simulations demonstrate that, beyond the left border of the stability 
region in Fig. 3(a) and beyond the lower border of the stability region in Fig. 
3(b), any initially created pulse decays to zero. This is quite understandable, 
as in the case when the gain coefficient 71 is too small, nothing but the trivial 
solution may exist as a steady state, i.e., no stationary pulse is present in this 
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Figure 3: Stability regions of the stationary two-dimensional solitary pulse in 
the parametric planes (71,72) for /3 = (a), and (/3,7i) for fixed 72 = 1/2 (b). 



case. On the other hand, various transitions to spatially extended patterns take 
place beyond the right border of the stability region in Fig. 3(a) and beyond the 
upper boundary in Fig. 3(b). For example, a persistent spatially extended and 
temporally chaotic state sets in for 71 > 2.6 at /3 = and 72 = 1/2. A transition 
to a localized zigzag (stepped) stripe state (discussed in the next section) takes 
place for 71 > 2.3 at /3 = -0.4 and 72 = 1/2. 

To check the effect of the third-order dispersion (TOD), a few numerical 
simulations of the accordingly modified equation ( pj]) were performed. Figure 5 
shows a snapshot of the 2D localized pulse affected by TOD with /3 = 0.1. The 
pulse is more circular than the one shown for the same values of the parameters, 
except that /3 = 0, in Fig. 1, and it is moving along the r direction. Figures 
6 (a) and (b) display the amplitude profiles a and the phase profiles (j) along 
two principal cross sections, x — L/2 and r — T,„ax, at the intersection of 
which the amplitude a{x,T) has its maximum (r^ax is not fixed as the pulse is 
moving). The width of the amplitude profile along the r direction is smaller than 
that for /3 = shown in Fig. 2(a), plausibly because the additional dispersion 
provided by TOD makes the pulse structure sharper. The phase profile along 
the X direction is strongly different from that in Fig. 2(b). Note also that the 
phase profile along the r direction is slightly asymmetric, which is a natural 
consequence of the symmetry breaking by the TOD term. 

3 One-dimensional localized pulses and the zigzag 
instability 

The 2D localized pulse is not the single stable pattern generated by Eq. (|l0|). 
Simulations also reveal a stable ID pattern in the form of a straight stripe 
localized in the cc-direction and uniform along the r axis. As it was mentioned 
in the previous section, this type of ID pulses can be found when the 2D pulse 
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^1 

Figure 4: The widths of the stable stationary two-dimensional pulse in the x- 
direction (rhombuses) and r-direction (crosses) vs. the cubic gain 71, at fixed 
values (i — Q and 72 ~ 1/2. 



becomes unstable (but it can also stably coexist with a 2D pulse, see below). 
In terms of the underlying optical-waveguide model, this solution corresponds 
to a spatial soliton, rather than a spatiotemporal one corresponding to the 
2D pulse considered in the previous section. Note that Eq. (|l0| ) is invariant 
against a Galilean-like transformation in the spatial domain with an arbitrary 
real parameter q, 

u{z, X, t) = exp ^iqx — 2^1^^ (-^^i 2;, r) , x = x ~ qz, (12) 

which implies a possibility of an arbitrary orientation of the spatial stripe soliton 
in the (z, x) plane. 

Obviously, there also exist more general steady-state solutions to Eq. ( p^ 
in the form of a stripe oriented at oblique directions in the spatiotemporal 
(x,r) plane, similar to the solution ^ to Eq. (|l|). In the optical waveguide, 
they correspond to moving spatial solitons. These oblique stripe patterns will 
appear below. 

Extended simulations demonstrate that the spatial solitons (stationary stripes) 
and 2D spatiotemporal pulses coexist, as stable solutions to Eq. ([lO[), in a 
broad parametric region, provided that the dispersion coefficient /? is positive 
(recall this corresponds to the case of anomalous chromatic dispersion in the 
waveguide). As /3 decreases, the stripes lose their stability against perturba- 
tions which rearrange them into a zigzag pattern. A typical example of a 
stable zigzag (which coexists with a stable 2D pulse) is shown in Fig. 7 for 
P = —0.05, 7i — 2.3 and 72 = 1/2. Both zig and zag segments of the pat- 
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Figure 5: A snapshot of a two-dimensional pulse obtained from Eq. (pi]), that 
includes the third-order-dispersion term with the coefficient f3 — 0.1, the other 
parameters being the same as in Fig. 1 . The absolute value of the field is shown 
as a function of x and r. 



tern are nothing else but pieces of the above-mentioned oblique stripes oriented 
under equal but opposite angles relative to the straight stripes. Note that the 
transformation (|l^) may as well be applied to the zigzag solution, generating 
zigzags with an arbitrary orientation in the (z, x) plane. 

While stable straight and zigzag stripes coexist with the stable 2D pulses, we 
have found that they never coexist with each other as stable solutions. Roughly, 
the straight stripes are stable at /3 > and unstable at /3 < 0. A bifurcation 
which destabilizes the straight stripe and gives rise to the zigzag is a supercritical 
(forward) one, see a typical example in Fig. 8. As a parameter characterizing 
the transition to the zigzag, in this figure we show the absolute value of the 
slope of the step's oblique segments (see Fig. 7) vs. /3 for L = 20. 

If the system's size Lr in the r direction is large enough, the instability of 
the straight stripe initiated by random perturbations may at first give rise to 
multi-zigzag patterns. However, a coarsening process occurs, and as a result 
the pattern eventually relaxes to a single-zigzag state, as is shown in Fig. 9 for 
Lj; = 20 and Lr = 60 at f3 = —0.05, 71 — 2.3 and 72 = 1/2. It means that 
two oblique stripes with the absolute value of the slope given by the bifurcation 
diagram shown in Fig. 8 are stable solutions in the case when the straight stripe 
is unstable. 

Zigzag patterns are known in models of the complex-Swift-Hohenberg po| 
and Segel-Newell- Whitehead (SNW) jljl types; coarsening of the zigzag struc- 
ture is also a known feature of the SNW equation. In particular, in the former 
model it was found that roll-type standing-wave patterns give rise to two spe- 
cific modes of small perturbations around them: one mode accounts for zigzag 
corrugation of the stripes, and the other one is a phase modulation of the os- 
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Figure 6: The amplitude (a) and phase (b) profiles of the pulse shown in Fig. 5, 
taken along the cross sections r = Tmax (continuous lines) and x = L/2 (dashed 
lines). The r axis is offset by L/2 — Tmax to stress that the peak positions of 
the amplitude and phase profiles coincide. 



dilations. The two modes are decoupled in the linear analysis, but they get 
coupled in a nonlinear regime. The same set of two modes can be found around 
the localized straight-stripe pattern in the present model, hence the onset of 
the zigzag instability in the present model is similar to that studied in detail in 
Refs. [Q. However, physical interpretation of each zigzag in terms of the pla- 
nar optical waveguide underlying the present model is quite different: a spatial 
soliton moving in the transverse direction in the homogeneous planar waveguide 
suddenly reverses its velocity to the opposite value. This appears to be the first 
example of such a "boomerang" motion mode of a stripe soliton in nonlinear 
optics (note the GL equation does not conserve momentum). 

4 Conclusion 

In this work, we have introduced a model of a two-dimensional optical waveguide 
with the Kerr nonlinearity, linear and quintic losses, cubic gain, and temporal- 
domain filtering. In the general case, the model also contains temporal disper- 
sion, although the latter ingredient is not necessary. The model may provide 
for a realistic average description of a nonlinear planar waveguide incorporated 
into a closed optical cavity. It takes the form of a two-dimensional cubic-quintic 
Ginzburg-Landau equation with an anisotropy of a novel type: the equation 
is diffractive in one direction, and diffusive (or mixed diffusive/dispersive) in 
the other. Following analogy with the model's one-dimensional counterpart, we 
have demonstrated by systematic direct simulations that the model gives rise 
to stable fully localized two-dimensional pulses, which, in terms of the optical 
planar waveguide, are spatiotemporal "light bullets" , existing due to the simul- 
taneous balances of diffraction, dispersion, and Kerr nonlinearity, and of linear 
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Figure 7: An example of a quasi-one-dimensional stable stationary pattern con- 
taining a zigzag. The absolute value of the field is shown as a function of x and 
r. This solution to Eq. (0) was found at /? = -0.05, 71 = 2.3, 72 = 1/2. 



losses (including those induced by the filtering), cubic gain, and quintic loss. 
A stability region of the two-dimensional pulses was identified in the system's 
parameter space. We have also found that the model generates stable (quasi-) 
one-dimensional patterns in the form of a simple straight stripe or oblique stripes 
with zigzags. The straight and oblique stripes may stably coexist with the two- 
dimensional pulse, but not with each other. A bifurcation which destabilizes the 
straight stripe and gives rise to the oblique ones was found to be supercritical. 
In terms of the optical waveguide, the straight-stripe solution corresponds to a 
simple spatial soliton, while the zigzag on oblique stripes may be realized as a 
spatial soliton moving in the transverse direction which suddenly changes the 
sign of its velocity. It seems to be the first example of such a dynamical behavior 
in nonlinear optics. 
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Figure 8: An example of the forward bifurcation that destabihzes the straight- 
stripe solution and gives rise to a stable zigzag pattern. In this case, 71 = 2.3 
and 72 — 1/2. The absolute value of the slope of the zigzag's oblique sections 
is plotted as a function of the dispersion coefficient (3. 
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